Preprint typeset in JHBP style - HYPER VERSION 



UH-511-1202-12 



Post-LHC? fine-tuning in the mSUGRA/CMSSM 
model with a 125 GeV Higgs boson 



Howard Baer", Vernon Barger^, Peisi Huang^, Dan Mickelson", Azar Mustafayev^ and 
Xerxes Tata"^ 

"^Dept. of Physics and Astronomy, University of Oklahoma, Norman, OK 73019, USA 

^Dept. of Physics, University of Wisconsin, Madison, WI 53706, USA 

"Dept. of Physics and Astronomy, University of Hawaii, Honolulu, HI 96822, USA 



E-mail: 


baerOnlm . ou . e&u 




bargerOpheno . wise . edu 


, IphuangTOwisc . edu 


dsmickelsonOou . edu 




azarOphys . hawaii . edi: 




tataOphys .hawaii . edu 



D. 

^ ■ Abstract: The recent discovery of a 125 GeV Higgs-like resonance at LHC, coupled with 

the lack of evidence for weak scale supersymmetry (SUSY), have severely constrained SUSY 
models such as mSUGRA/CMSSM. As LHC probes deeper into SUSY model parameter 
space, the little hierarchy problem - how to reconcile the Z and Higgs boson mass scale with 
the scale of SUSY breaking ~ will become increasingly exacerbated unless a sparticle signal 
is found. We evaluate two different measures of fine-tuning in the mSUGRA/CMSSM 
model. The more stringent of these, AhS) includes effects that arise from the high scale 
origin of the mSUGRA parameters while the second measure, Aew? is determined only by 
weak scale parameters: hence, it is universal to any model with the same particle spectrum 
and couplings. Our results incorporate the latest constraints from LHC7 sparticle searches 
and require a light Higgs scalar with rrih ~ 123 — 127 GeV. We present fine-tuning contours 
in the mo vs. mii2 plane for several sets of Aq and tan/3 values. We also present results 
for Ahs and Aew from a scan over the entire viable model parameter space. We find a 
Ahs ^ 10^, or at best 0.1% fine-tuning. For the less stringent electroweak fine tuning, 
we find Aew ^ 10^, or at best 1% fine-tuning. Two benchmark points are presented 
that have the lowest values of Ahs and Aew- Our results provide a quantitative measure 
for ascertaining whether or not the remaining mSUGRA/CMSSM model parameter space 
is excessively fine-tuned, and so could provide impetus for considering alternative SUSY 
models. 
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1. Introduction 



The recent spectacular runs of LHC at ^/s = 7 and 8 TeV have led to identification of 
a Higgs-like boson^ with mass rrih ~ 125 GeV [j^, ^. This is in accord with predictions 
from the minimal supersymmetric standard model (MSSM) which requires that the lighter 
higgs scalar mass m/j < 130 — 135 GeV Q. Since values of mh > Mz are only possible 
due to radiative corrections, the upper end of the range depends on the masses of third 
generation sparticles that one is willing to allow. To achieve rrih ~ 125 GeV, either large 
mixing or several TeV masses are required in the top squark sector. In models such as 
the much-studied minimal supergravity (mSUGRA or CMSSM) model [Q, |5|, values of 
trilinear soft breaking parameter |^o| ~ (1-5 — 2)mo are favored, along with top squark 
masses m^^ ^>1 — 2 TeV: for positive values mo is typically larger than 5 TeV ^ 

While the measured value of nih is within the expected range of even the simplest SUSY 
models, there is at present no sign of SUSY particles at LHC. From LHC data analyses 
within the mSUGRA model, mass limits of > 1.4 TeV when rriq ~ nig and rrig > 0.9 TeV 
when rriq S> rrig have been reported[Q, |9[. Several groups [10| have updated their fits of the 
mSUGRA/CMSSM model to various data sets, now including information from LHC7 and 
LHC8 Higgs-like boson discovery and LHC7 sparticle mass limits. Typically, the best fit 
regions have moved out to large values of mo and mi/2 to accomodate the LHC sparticle 
mass limits and Higgs discovery. Such large mo and mi/2 values lead to sparticle masses in 
the multi-TeV mass range, thus exacerbating what has become known as the little hierarchy 
problem: how do such large SUSY particle masses and soft breaking parameters conspire 
to yield the weak scale typified by the Z-boson mass Mz — 91.2 GeV. The conflict between 
the strong new LHC sparticle mass limits and the comparatively low values of Mz and nih 
has intensified interest in the fine-tuning in supersymmetric models|ll, 12, |l^, 



To set the stage for this analysis, we begin by reviewing radiative corrections (assumed 
perturbative) to scalar field masses. In a generic quantum field theory, taken to be the 
low energy effective theory whose domain of validity extends up to the energy scale A, the 
physical mass squared of scalar fields takes the schematic form (at leading order), 

,2 _ ^2 I 9^ a2 I 9^ ^2 f \ I 5^ ^2 

how 



^4> = '^'PO + ^lYfo?^ + ^2Y^"^iowlog ( 1 + Csj^m,^^. (1.1) 



In Eq. (1.1), g denotes the typical coupling of the scalar cf), m^g is the corresponding mass 
parameter in the Lagrangian, IGvr^ is a loop factor, and Ci are constants that aside from 
spin, colour and other multiplicity factors are numbers 0{1). The scales mio„ and A re- 
spectively denote the highest mass scale in the effective theory and the scale at which this 
effective theory description becomes invalid because heavy degrees of freedom not included 
in the low energy Lagrangian become important. For instance, if we are considering cor- 
rections to the Higgs sector of the MSSM is embedded into a Grand Unified Theory (GUT) 
framework, A ~ Mqut and miow ~ -^SUSY (or more precisely, miow is around the mass of 



^This particle has spin or > 2 and couples directly to the ZZ, and with weaker evidence also to the 
WW, systems. The latter property implies a connection with electroweak symmetry breaking, characteristic 
of the Higgs boson. 
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the heaviest sparticles that have large couphngs to the scalar (j)). Finally, the last term in 
comes from loops of particles of the low energy theory, and their scale is set by miow 
These terms may contain logarithms, but no large logarithms since effects of very high 
momentum loops are included in the Ci and C2 terms. These finite corrections provide 
contributions to that which we have referred to as electroweak fine-tuning in a previous 
study PI. 



If the effective theory description is assumed to be valid to the GUT scale, the Ci 

2 

</>o 



term is enormous. Even so it is always possible to adjust the Lagrangian parameter w? 



to get the desired value of m^2 < m^^^- This is the hig fine tuning problem of generic 
quantum field theory with elementary scalars. This problem is absent in softly broken 



supersymmetric theories because Ci = 0. We see from Eq. (IT) that if the physical value 
of nrifi) is significantly smaller than miow (which in the case of the MSSM ~ "^tj); we will 
still need to have significant cancellations among the various terms to get the desired value 
of m^. This is the little hierarchy problem. We also see that in models such as mSUGRA 
that are assumed to be valid up to very high energy scales A ~ Mqut — Mp, the magnitude 
of the C2 term typically far exceeds that of the C3 term because the logarithm is large, 
and hence is potentially the largest source of fine-tuning in such SUSY scenarios. 

Because the C2 and C3 terms in Eq. ( |1.1| ) have somewhat different origins - the C2 
term represents corrections from physics at scales between miow and A, while the C3 term 
captures the corrections from physics at or below the scale miow ~ we will keep individual 
track of these terms. In the following we will refer to fine-tuning from C2 type terms as 
high scale fine-tuning (HSFT) (since this exists only in models that are valid to energy 
scales much larger than miow) and to the fine tuning from Cs-type terms as electroweak 
fine-tuning (EWFT) for reasons that are evident. We emphasize that the sharp distinction 
between these terms exists only in models such as mSUGRA that are assumed to be 
a valid description to very high scales, and is absent in low scale models such as the 
phenomenological MSSM [16|. 

In this paper, we quantify the severity of fine-tuning in the mSUGRA model, keeping 
separate the contributions from the two different terms. We are motivated to do so for two 
different reasons. 

• First, as emphasized, C2 type terms appear only if the theory is applicable out to 
scale A ^ miow 7 while the C3 type terms are always present. In this sense, the fine- 
tuning from the C3 type terms is ubiquitous to all models, whereas the fine-tuning 
associated with the (potentially larger) C2 type terms may be absent, depending on 
the model. 

• Second, as we will explain below, there are two very different attitudes that one can 
adopt for the fine-tuning from C2 type terms. Keeping the contributions from C2 
and C3 separate will allow the reader the choice as to how to interpret our results 
and facilitate connection with previous studies. 

The remainder of this paper is organized as follows. In Sec. ^ we introduce our measures 
of fine-tuning. As usual, we adopt the degree to which various contributions from the 
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minimization of the one-loop effective potential in the MSSM Higgs boson sector must 
cancel to reproduce the observed value of M| as our measure of fine tuning. We use these 
considerations to introduce two different measures. The first of these is the less stringent 
one and relies only on the weak scale Lagrangian that arises from mSUGRA with total 
disregard for its high scale origin, and is referred to as electroweak fine-tuning (EWFT). The 
other measure that we introduce incorporates the high scale origin of mSUGRA parameters 
and is therefore referred to as high scale fine-tuning (HSFT). In Sec. ^ we present contours 
for both HSFT and EWFT in several mSUGRA mo vs. mi/2 planes along with excluded 
regions from LHC7 sparticle searches. We find that while LHC7 sparticle mass limits 
typically require EWFT at ~ 1% level, the requirement that ruh ~ 125 GeV requires much 
more severe EWFT in the 0.1% range in the bulk of parameter space. As anticipated, 
HSFT is even more severe. We also find that the hyperbolic branch/focus point region 
(HB/FP)p7|] - while enjoying lower EWFT than the bulk of mSUGRA parameter space - 
still requires fine-tuning at about the percent level. The fine-tuning situation is exacerbated 
by the requirement of large IAq/'^oI for which the HB/FP region is absent, resulting in 
large EWFT (and even larger HSFT). In Sec. ^ we present results from a complete scan 
over mSUGRA/CMSSM parameter space. In this case, respecting both the LHC7 sparticle 
mass bounds and ruh = 123 — 127 GeV (in accord with the estimated theory error on our 
calculation of nih), we find parameter space points with maximally 0.1% HSFT and 1% 
EWFT. We leave it to the reader to assess how much fine-tuning is too much, and also to 
judge the role of HSFT in models such as mSUGRA/CMSSM that originate in high scale 
physics. We present and qualitatively discuss the phenomenology of two model points with 
the lowest HSFT and the lowest EWFT in Sec. |. We end with some concluding remarks 
and our perspective in Sec. ^ 

2. Fine-tuning 

We begin by first writing the Higgs potential whose minimization determines the elec- 
troweak gauge boson masses as, 

Vmggs = {ml^ + f^')K\^ + + f^')\hT 

-B^il^A + h.c.) +^-{g' + g"){\hlf - \h'ff + AV , (2.1) 

where the radiative corrections (in the one-loop effective potential approximation) are given 
in the DR scheme by, 

Here, the sum over i runs over all fields that couple to Higgs fields, is the Higgs 

field dependent mass squared matrix (defined as the second derivative of the tree level 
potential), and the trace is over the internal as well as any spin indices. One may compute 
the gauge boson masses in terms of the Higgs field vacuum expectation values Vu and Vd by 



log 



MjM 



(2.2) 
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minimizing the scalar potential in the and directions. This leads to the well-known 
condition 

ml _ «, + S^)-(rnk+S^)tan2/3 , 

~" tan2/3-l ^ • ^^-^^ 

Here the and Ti'^ terms arise from first derivatives of /SV evaluated at the potential 
minimum and tan/3 = Vu/vd- At the one-loop level, contains the contributions £^(^1^2)5 
S:i(6i,2), S^(fi,2), S^(H?i,2), S^(^i-4), and 

contains similar terms along with and SJ^(t) while = [ p^ . 

Although we have highlighted third generation matter sfermion contributions here 
because these frequently dominate on account of their large Yukawa couplings, we note 
that there are also first/second generation contributions Ti^{q^t) and that arise 

from the quartic D-term interactions between the Higgs sector and matter scalar sector 
even when the corresponding Yukawa couplings are negligibly small. These contributions 
are proportional to (T3. — Qjsin^^^y) x F(m?), where T3. is the hypercharge, Qj is the 
electric charge and F{m?) = m^(log ^ — 1) of the i^^ matter scalar. Although the scale of 
these is set by the electroweak gauge couplings rather than the top Yukawa coupling, these 
can nevertheless be sizeable if the squarks of the first two generations are significantly 
heavier than third generation squarks. However, in models such as mSUGRA - where 
all squarks of the first two generations (and separately, the corresponding sleptons) are 
nearly mass degenerate - these contributions largely cancel. Indeed, the near cancellation 
(which would be perfect cancellation in the case of exact degeneracy) occurs within each 
generation, and separately for squarks and for sleptons. These terms, summed over each of 
the first two generations, are always smaller than the other terms in the Cj and Bi arrays 
used to define our fine-tuning criterion below, and so do not alter our fine-tuning measure 
defined below. 

The reader may wonder that we are treating the first two generations differently from 
the third generation in that for the latter we consider the contributions from each squark 
separately (i.e. not allow for cancellations of the contributions to say from different 
squarks to cancel), while we sum the contributions from the entire first/second generation 
to obtain a tiny contribution. The reason for this is that the mSUGRA framework predicts 
degenerate first/second generation squarks (and sleptons) while the top squark masses 
(remember that top squarks frequently make the largest contribution to S") are essentially 
independent. In an unconstrained framework such as the pMSSM we would not combine 
the contributions from the first /second generation scalars; if these are very heavy and have 
large intra-generation splitting, their contribution to Aew can be significant. 



2.1 Electroweak scale fine-tuning 



One measure of fine-tuning, introduced previously in Ref. [14, 12], is to posit that there are 
no large cancellations in Eq. (|2.3|) . This implies that all terms on the right-hand side to be 



comparable to M|/2, i.e. that each of the three tree level terms Ch^ = l"^lf^/(tan^ /3 — 1)|, 
= I ~ tan^ /3/(tan^ /3 — 1)|, = | — fi'^\ and each C^u.d... is less than some 

characteristic value A where A ~ M^. (Here, i labels SM and supersymmetric particles 
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that contribute to the one-loop Higgs potential and includes the sum over matter sfermions 
from the first two generations.) This leads to a fine-tuning measure 



max{Ci)/{Ml/2) 



(2.4) 



A feature of defining the fine-tuning parameter solely in terms of weak scale parameters 
is that it is independent of whether the SUSY particle spectrum is generated using some 
high scale theory or generated at or near the weak scale, as in the pMSSM or possibly in 
gauge- mediation |]l8|: if the spectra and weak scale couplings from two different high scale 
theories are identical, the corresponding fine-tuning measures are the same. However, as we 
will see in Subsection |2[^, in theories such as mSUGRA Aew does not capture the entire 
fine-tuning because Eq. ( p.3| ) does not include information about the underlying origin of 
the weak scale mass parameters. 

It is worthwhile to note that over most of parameter space the dominant contribution 
to Aew comes from the weak scale values of Tn?H^ and . To see this, we note that 
unless tan /3 is very small, aside from radiative corrections, we would have simply that 



M|/2 



-m 



/i^. As is customary, the value of /i^ is selected so that the correct value 



of is generated. In this case, over much of parameter space Aew ~ 1^^1/(^1/2)- Only 
when 1^1 becomes small do the radiative corrections become important - providing the 
largest contribution to Eq. (|2.3|). Thus, contours of fixed Aew typically track the contours 
of |//| except when |//| is small; in this latter case, Aew is determined by the whose 
value is loop-suppressed. In Fig. |l| we show the surface of values in the mg vs. mi/2 
plane of mSUGRA/CMSSM for = and tan/3 = 10. Here, /i is small either at low ttiq 
and mij2 (the bulk region [p^]), or in the HB/FP region[17| at large values of mo- 




m,/2 [TeV] 



Figure 1: The value of fi in the mo vs. mi/2 plane of mSUGRA for = and tan/3 = 10. We 
set /i = in theoretically forbidden regions. 
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2.2 High scale fine-tuning 

As mentioned above, Eq. (|2.3| ) is obtained from the weak scale MSSM potential and so 
contains no information about its possible high scale origin. To access this, and make 
explicit the dependence on the high scale A, we must write the weak scale parameters 



mjj- in Eq. (2^) as 



,2 



m 



where m?^ ^{A) and /i^(A) are the corresponding parameters renormalized at the high scale 
A. It is the 5m|^ terms that contain the log A dependence shown in the C2 type terms 



in Eq. (IT). In this way, we get 



ml (m% (A) + dm% + S^J) - (A) + 6m% + S^J) tan^ /3 

^ = ^ ' ^ ^ — \ ' ^ --{|?{K) + 5^l^) . (2.5) 

2 tan p — 1 

Following the same spirit that we had used in our earlier analyses [^], we can now define 
a fine-tuning measure that encodes the information about the high scale origin of the 



parameters by requiring that each of the terms on the right-hand-side of Eq. (2^) to be 
smaller than a pre-assigned Ahs times The high scale fine-tuning measure Ahs is 

thus defined to be 

Ahs = max{Bi)/{Ml/2) , (2.6) 

with 

Bh, = |m2,^(A)/(tan2/3 - 1)|, Bsh, = \5mlj{t&n^^- 1)|, 
Bh^ = I -ml^JA)tan2/3/(tan2/3- 1)1, Bsh^ = | - ^m^^^ tan^ /^/(tan^ /? - 1)|, etc., 



defined analogously to the set Cj in Sec. 2.1. As discussed above, in models such as 



mSUGRA whose domain of validity extends to very high scales, because of the large log- 
arithms one would expect that (barring seemingly accidental cancellations) the B^u^ con- 
tributions to Ahs would be much larger than any contributions to Aew because the rn^H^ 
evolves from mg to negative values. 

As we have noted, Aew indeed provides a measure of EWFT that is determined only 
by the sparticle spectrum: by construction, it has no information about any tuning that 
may be necessary in order to generate a given weak scale SUSY mass spectrum. Thus, 
while for a given SUSY spectrum Aew includes information about the minimal amount 
of fine-tuning that is present in the model, Ahs better represents the fine-tuning that is 
present in high scale models. 

The reader may have noticed that - unlike in our definition of Aew in Eq. (^3) where 
we have separated out the contributions from various sources and required each of these 
to not exceed some preassigned value - we have neglected to separate out the various 
contributions to ^ that determine Ahs- We have done so mainly for convenience,^ 

but this will also help us to connect up with what has been done in the literature. 



^Unlike for Aew where we have separated the contributions by particles (and treated these as indepen- 
dent) for the electroweak scale theory, in a constrained high scale model, these would not be independent. 
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Before closing this section, we remark that our definition of Ahs differs in spirit from 
that used by some groups [|l5|. These authors write the rnjj^ as a quadratic function of the 
high scale parameters = {mo, mi/2, Aq} for mSUGRA, i.e. 

^Hn = "-ij^i^j ' (2-7) 

and substitute this (along with the corresponding form for m'jj^) in Eq. ( |2.3D to examine 
the sensitivity of to changes in the high scale parameters.^ In the resulting expression, 
the coefficient of ttiq in Eq. ( |2.7D is often very small because of cancellations with the 
large logarithms, suggesting that the region of mSUGRA with rather large ttt-q (but small 
mi/2 ^-iid Aq) is not fine-tuned: we feel that this is misleading and so have separated 
the contributions from the large logarithms in our definition of Ahs- Combining all rriQ 
contributions into a single term effectively combines m?(A) + 5mf into a single quantity 
which (aside from the one-loop terms and S^) evidently is the weak scale value of m? 
in our definition of Ahs- Except for these one-loop correction terms, Ahs then reduces to 
Aew! 

In defining Ahs as above, we have taken the view that the high scale parameters as 
well as the scale at which we assume the effective theory to be valid are independent. In the 

absence of an underlying theory of the origin of these parameters, we regard cancellations 



between terms in Eq. (2^) that occur for ad hoc relations^ between model parameters and 
lead one to conclude that Mz is not fine-tuned as fortuitous, and do not incorporate it 
into our definition of high scale fine-tuning. We emphasize that we would view the fine- 
tuning question very differently if indeed the high scale parameters were all related from 
an underlying meta-theory.^ In that case, though, as we just mentioned, Aew would be 
an adequate measure of fine tuning. 



3. Results in rriQ vs. mi/2 plane 

We present our first results as contours of Ahs and Aew ™ the tuq vs. plane of 

the mSUGRA/CMSSM model. For all plots, we take m-j = 173.2 GeV and we generate 
SUSY particle mass spectra from Isasugra v7.83 [^. In Fig. ||, we show contours of Ahs 
in frame a) and for Aew in frame b). For both frames, we take = 0, tan/3 = 10 and 
> 0. The gray-shaded regions running from the extreme left of the plot, across the 
bottom and on to the right are excluded by either a fi as LSP (left-side), LEPl constraints 

Instead, we could separate out contributions that have independent origins in the high scale model. For 
instance, for the mSUGRA model we should separately require contributions from gauginos, scalars and 
A-parameters to 5m%. to be small. We have not done so here mainly for expediency. In this sense if 
accidental cancellations reduce Ahs to very small values, this should be interpreted with care. 

^Typically these authors use A = -f^-^^ (where labels the input parameters) as a measure of 
the sensitivity to parameters This prescription agrees with our A at tree level, but differs when loop 
corrections are included. 

*It may be argued that such an analysis is helpful as a guide to model builders attempting to construct 
models of natural SUSY. 

^This situation seems to occur in the so-called mixed- modulus-anomaly mediated SUSY breaking models 
for some ranges of the mixing parameter a as emphasized in Ref. ]2c|]. 
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(bottom) or lack of appropriate EWSB (right-side). The region marked LEP2 is excluded 
by LEP2 chargino seaches {tn^^^ > 103.5 GeV) |2^. The region below the contour labeled 
LHC7 is excluded by lack of a SUSY signal from SUSY searches at LHC7 with 5 fb~^ 
of data|8|, ^. The dashed portion of the contour is our extrapolation of LHC7 results 
to higher values of mo than are shown by the Atlas/CMS collaborations. The green- 
shaded region is where the thermally-generated relic density of neutralinos (computed 
usmg IsaReD[||]) satisfies VL^^y < 0.1194, the 2a upper limit on the density of cold dark 



matter obtained by the WMAP collaboration |24|. This region encompasses the stau- 



coannihilation strip |^5[ (extreme left), the bulk region ]19[] (bottom left corner) and the 



well-known focus point/hyperbolic branch region [17] of the model. The shaded region 



labeled is where the measured muon magnetic moment [^] satisfies 4.7x10 < < 

-10 



52.7 X 10~ , within 3a of its theoretical value [27|. For Aq = adopted in this figure, 
mil < 123 GeV over the entire parameter plane, so that mSUGRA is excluded for ~ 
(as noted in Ref. p) unless one has very high values of mo and mi/2 p^]. 

As might be anticipated, Ahs grows with increasing values of mo or mi/2, so that we 
expect contours of fixed Ahs to be oval-shaped in the mo — mi/2 plane. This is readily seen 
in frame o) of Fig. ^, except that because the oval is extremely elongated since the scales 
on the two axes are very different, we see only a small part of this contour (which appears 
as nearly vertical lines) for very large values of Ahs- We have checked that Ahs < 150 is 
already excluded by LHC searches, so high scale fine-tuning of less than a percent is now 
mandatory for ^o = 0. If we take the high scale origin of the mSUGRA model seriously, we 
see that without a theory that posits special relations between the parameters that could 
lead to automatic cancellation of the large logarithms that enter Ahsj we are forced to 
conclude that LHC data imply that the theory is fine-tuned to a fraction of a percent. For 
the portion of the plane compatible with LHC constraints on sparticles, the smallest values 
of Ahs occur where mo and mi/2 are simultaneously small. As mo moves to the multi-TeV 
scale, Ahs exceeds 1000, and fine-tuning of more than part per mille is required. 

In frame b) of the figure, we show contours of constant Aew- Over most of the plane, 
these contours tend to track contours of constant ^'^ since M^/2 ~ —m'jj^—^'^ so that when 
|m|^^| ^ M^/2, then —m'jj^ ~ ^u^. Thus, along the contours of Aew? the value of m^^ is 
independent of mg at least until the contours turn around at large values of mg and mi/2- 
This is just the focus point behaviour discussed in the second paper of Ref. p^.^ The Aew 
contours, for large values of mo bend over and track excluded region on the right where fj? 



becomes negative. This is the celebrated hyperbolic branch |17] of small The contours 
of Aew then bend around for very large values of mo because contributions, especially 
from ^2 loops — increase with mo — begin to exceed —m?^^ ~ /U^. Indeed, Fig. ^5) shows 
that there is a region close to (but somewhat removed from) the "no EWSB" region on the 
right where Aew becomes anomalously small even for large values of mo and mi/2- It is 
instructive to see that while this low EWFT region is close to the relic-density consistent 



®More precisely, the discussion in this paper was for a fixed value of mi/2 so that the range of mo was 
limited because we hit the theoretically excluded region. We see though that the same value of can be 
obtained if we simultaneously increase mo and mi/2 so that we remain in the theoretically allowed region. 
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region with small /x [|^, it is still separated from it J While Aew ~ 100 is excluded at 
low mo, this 1% EWFT contour, even with the resolution of our scan, extends out to very 
large uiq ~ 6 TeV values for mi/2 as high as 1 TeV! While these plots show that relatively 
low EWFT (Aew of a few tens) is still allowed by LHC7 constraints on sparticles, it is 
important to realize that these planes are now excluded since they cannot accommodate 
nih ~ 125 GeV. 



mSUGRA: tanf,=10, 4„ =0, n >0, m, =173.2 GeV mSUGRA: tanf,=10, A„ =0, \i >0, m, =173.2 GeV 




0123456789 10 0123456789 10 



mo (TeV) m,(TeV) 

Figure 2: Contours of a) Ahs and b) Aew in the mSUGRA model with Aq — and tan/3 — 10. 
We take ^ > and mt = 173.2 GeV. The grey region on the left is excluded either because the stau 
is too light or becomes tachyonic, the grey region at the bottom is excluded by LEPl constraints, 
while in the region on the right we do not get the correct pattern of EWSB, since either /i^ or 
m\ become negative. The region labeled LEP2 is excluded by constraints on the chargino mass. 
The region labeled is allowed at the 3a level by the E821 experiment while in the green-shaded 
region, the thermal neutralino relic density is at or below the WMAP measurement of the cold dark 
matter density. The lighter Higgs boson mass < 123 GeV throughout this parameter plane. 

Before moving on to other planes, we remark that for the smallest values of mo in the 
LHC-allowed regions of the figure, Ahs ~ Aew- As we have explained, Ahs is determined 
by the value of |(5m|^^| (see Eq. ^), which for mo ~ is just I^t-^^^I that determined Aew 
when mo is very small. We thus see that the two measures are roughly comparable for small 
values of mo but deviate from one another as mo is increased. We see that Ahs typically 
exceeds Aew by an order of magnitude, because of the large logarithm of the ratio of the 
GUT and weak scales, except in the HB/FP region where Aew is exceptionally small. 

In Fig. ^ we show the mo vs. mij2 plane for tan/3 = 50 and Aq = 0. The contours 
in both frames are qualitatively very similar those for the tan /3 = 10 case. As expected, 
regions of low Aew extend to very large mo and mx/2 in the HB region. One difference 
from the tan /3 = 10 case discussed above is that this time the HB region largely overlaps 

^Much of the literature treats these regions as one. While this is fine for some purposes, it seems 
necessary to be clear on the difference when discussing either dark matter or EWFT. Note that Ahs is 
large in both regions. 
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with the rehc-density-consistent green-shaded region. Again, the entire region of plane 
shown is excluded by the LHC Higgs discovery at 125 GeV. 



mSUGRA: tanf,=50, 4„ =0, n >0, m, =173.2 GeV 
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Figure 3: Contours of a) Ahs and b) Aew in the mSUGRA model with Aq — 0, tan /3 — 50 and 
/i > 0. The value of mt as well as the various shaded/coloured regions are as in Fig. |^. 

In Fig. we show contours of Ahs find Aew for tan/3 = 10 and Aq = —mo. The first 
thing to notice is that the HB/FP region does not appear. The region at extremely large 
niQ is still theoretically excluded, but more typically because m\ turns negative (or there 
are tachyons) not because /i^ turns negative.^ In addition, the very large tuq ^7 — 9 TeV 
region yields a value of ruh > 123 GeV: thus, the bulk of this plane is still excluded. The 
contours of Ahs are qualitatively similar to the = cases, and LHC7 still excludes 
Ahs < 100, so again a HSFT of more than 1% is required. In the region with nih > 
123 GeV, Ahs ^ 1-5 x 10^, and extreme HSFT is required. Moving to frame b), we note 
that though the contours of fixed Aew now run from top left to lower right, these still 
follow the lines of fixed values of /i^. Moreover, values of Aew below 100 are excluded 
by just the LHC7 sparticle mass constraints. If one also imposes rufi > 123 GeV, then 
Aew ^ 2000 is required over the entire plane shown. 

In Fig. ^, we show the rriQ vs. plane for Aq = —niQ but with tan/3 = 50. We see 
this is qualitatively very similar to the previous figure. Again the theoretically excluded 
region occurs at values of ttiq far beyond the range shown. Note though that the contour 
of ruh = 123 GeV has moved to slighly lower mo values. Still, requiring ruh > 123 GeV 
requires Ahs > 5 X 10^, and Aew > 700. 

According to Ref. Q, large mixing in the top squark sector and consequently the 
largest values of rn/i occur in mSUGRA for Aq ~ —2mQ. In Fig. ^, we show contours of 

*For mi/2 ~ 500 GeV, this happens for mo > 22 TeV. We mention that this breakdown of parameter 
space could be an artifact of the ISAJET algorithm for computing the sparticle mass spectrum in mSUGRA. 
An approximate tree-level spectrum is first required in order to evaluate the radiative corrections that can 
potentially yield a valid solution using an iterative procedure. But in the absence of a non-tachyonic, tree- 
level spectrum with the correct EWSB pattern, the program is unable to compute the radiatively corrected 
mass spectrum. 
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mSUGRA: fanp=JO, A„ =-m^ \i >0, m, =173.2 GeV 



mSUGRA: tan^=10, A„ =-m„, |i >0, m, =»73.2 Gel' 
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Figure 4: Contours of a) Ahs and b) Aew in the mSUGRA model with — ^fno, tan/? = 10 
and /i > 0. The value of mt as well as the various shaded/coloured regions are as in Fig. |[ 
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Figure 5: Contours of a) Ahs and h) Aew in the mSUGRA model with — ^fno, tan/? = 50 
and fi > Q. The value of mt as well as the various shaded/coloured regions are as in Fig. |[ 



Ahs and Aew for tan/3 = 10 and Aq = —2mQ. We note again that the HB/FP region 
does not appear in this plane. Notice also that the contours of m/j = 123 GeV have 
moved all the way down to mg ~ 2 TeV: thus, now much of the mSUGRA plane shown is 
allowed by the LHC Higgs-like resonance discovery. In fact, the portion of the plane with 
"^0 ^6 — 8 TeV gives too large a value of rrih > 127 GeV. The portion of the mo vs. mi/2 
plane allowed by both LHC7 sparticle searches and by having nih ~ 123 — 127 GeV requires 
Ahs ~ 10^ - 10^ or 0.1-0.01% HSFT. The EWFT required is Aew > 10^ also large. The 
lesson learned here is that the remaining mSUGRA regions with rrih ~ 123 — 127 GeV, and 
which obey sparticle mass constraints, are highly fine-tuned, even with the less restrictive 
EWFT measure. 
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mSUGRA: tanf,=10, Ag =-2m^ n >0, m, =173.2 GeV 



mSUGRA: tanf,=10, 4„ =-2m^ ^ >0, m, =173.2 GeV 




Fig. 1^ shows the mSUGRA plane for = — 2mo but with tan/3 = 50. In this case, 
large theoretically excluded parameter regions appear and these only grow larger until the 
entire parameter space collapses for even higher tan/3 ~ 55 — 60 12^. The region on the 
right is forbidden because m\ turns negative, not because |^| becomes small: this is why 
there is no DM-allowed region for large values of mo- The ruh = 123 GeV contour nearly 
coincides with Ahs = 10^ and Aew = 500. In this case, values of mo > 6 TeV are excluded 
as giving rise to too heavy a value of m^,. Thus, again the regions with nih ~ 123—127 GeV 
and obeying LHC7 sparticle search constraints, are highly fine-tuned. 
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Figure 7: Contours of a) Ahs and h) Aew in the mSUGRA model with — —2mo, tan/3 = 50 
and /i > 0. The value of nit as well as the various shaded/coloured regions are as in Fig. ^ 

Before closing this section, we digress to compare our results for the EWFT measure 
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with some results in the recent literature for the fine-tuning within the mSUGRA/CMSSM 
model calculated using the procedure described at the end of Sec. 2^. We have already 
argued at the end of that section that the fine-tuning measure that results from substitut- 
ing rrijj^ in Eq. ( |2.7| ) and the analogous expression for rrij^^ into Eq. ( |2.3D should match 
our EWFT measure. To check this, we have compared our results in Fig. ^b) to those in 
Fig. 1 of the first paper of Ref. [^]. There, these authors show the minimum value of 
their fine-tuning parameter A in the mo — Tni/2 plane, marginalizing over a range of 
and tan (3. We see that the shapes of their A contours are qualitatively similar (except in 
the large mo region where the contours turn around because radiative correction effects 
are important) to those of the contours in frame b) of Figs. |2| and ^ We use our Aq = 
figures for this comparison because of all the figures these have the smallest value of Aew- 
We have also checked that for any chosen value of mo and mi/2 ^ of Antusch et al. has a 
magnitude similar to (but never larger than) the corresponding lowest Aew that we obtain 
for any choice of Aq and tan /3. 



4. Scan over mSUGRA parameter space 

While the results of the previous section provide an overview of both the EWFT and 
the HSFT measures in light of LHC 7 and LHC8 constraints on sparticle and Higgs boson 
masses, we only presented results for particular choices of ^o and tan /3, and for /i > 0. In 
this Section, we present results from a scan over the complete mSUGRA parameter space 
with the following range of model parameters: 

mo : 0-15 TeV, (4.1) 

mi/2 : - 2 TeV, (4.2) 

-2.5 < Ao/mo{3) < 2.5, (4.3) 

tan/3 : 3-60. (4.4) 

We will show results for both fj, > and /U < 0. For each solution generated, we require 

1. electroweak symmetry be radiatively broken (REWSB), 

2. the neutralino Zi is the lightest MSSM particle, 

3. the light chargino mass obeys the LEP2 limit that > 103.5 GeV p^ , 

4. m/i = 125ib2 GeV in accord with the recent Higgs-like resonance discovery at LHC |l], 
^ and 

5. the mass spectra obey LHC7 sparticle mass constraints for the mSUGRA model |^. 

Our first results are shown in Fig. ^ for a) Ahs and b) Aew versus mg. Solutions with 
fi < are shown as blue circles while solutions with fi > are shown in red crosses. In frame 
a), we see that Ahs values occupy a rather narrow band which increases monotonically 
with mo. Values of mo < 1 TeV are excluded by the requirement ruh > 123 GeV. The fi > 
and /i < solutions occupy essentially the same band. This is not surprising because the 
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large logarithms are essentially independent of the sign of ^. The minimum allowed value 
of Ahs is ~ 1000, so that at least 0.1% fine-tuning is required of all remaining mSUGRA 
solutions. The minimum for Ajjs occurs at mo ~ 1500 GeV. This minimal Ahs solution is 
shown as a benchmark point in Sec. ^. For niQ as high as 15 TeV, Ahs increases to nearly 
10^. In frame 6), we show Aew versus ijiq. Here, the shape of the allowed region is very 
different from the Ahs case in frame a). Low values of can give nih > 123 GeV only 
if |j4o/mo| is sizeable and, as we have already seen, yield Aew of at least several hundred. 
Smaller values of Aew are obtained only in the HB/FP region where mg is large. In 
other words, in the "hole region" in frame b), we have mh < 123 GeV. The point with the 
minimum value of Aew ~ 100 occurs at mo ~ 7900 GeV, and is shown as the electroweak 
benchmark point in the next section. Over the remaining mSUGRA parameter space, at 
best 1% EWFT is required. 




Figure 8: Fine-tuning measures Ahs and Aew versus niQ from a scan over mSUGRA/CMSSM 
parameter space for /i > (red crosses) and /i < (blue circles). We take mt — 173.2 GeV. 

In Fig. ^, we show the distributions of Ahs and Aew versus mi/2- The sharp edge 
on the left is a reflection of the lower limit on rrig from LIIC7 searches. In frame a) for 
AhSi we see that the minimal Ahs is spread across a wide spectrum of mi/2 values. This 
is consistent with the behavior of Ahs shown in Fig. §a) where the HSFT contours are 
nearly vertical, indicating little dependence on mi/2- In frame b), the minimal values of 
Aew are also spread across the mi/2 range. For both Ahs and Aew, there may be a 
slight preference for lower mi/2 values. 

In Fig. |l^ we show how Ahs and Aew are distributed versus Aq/itiq. In frame a), 
we see that minimal Ahs ~ 1000 is obtained for Aq/uiq ~ —2, which is also the vicinity 
of where rrih is maximal for given mg and mi/2 values. There is also a minimum at 
AQ/rriQ ~ 2.5, with Ahs reaching only to ~ 3000 ^. In frame 6), the value of Aew is even 
more correlated with AQ/rriQ. For l^lo/mol ^ 1, Aew tends to be smaller than for larger 
values of |Ao/mo|. The solutions with the least EWFT occur at Ao/mo ~ ±0.6, with the 
minimal Aew ~ 100. Once again, this occurs in the HB/FP region mentioned above. For 
larger magnitudes of ^o/mo, the HB/FP region is absent, and Aew is much larger. The 

^The asymmetry of the minimum of Ahs with respect to the sign of Ao may only be a reflection of the 
fact that it is more difficult to generate large values of ruh for positive values of Aq. 
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Figure 9: Fine-tuning measures Ahs and Aew versus mi/2 from a scan over mSUGRA/CMSSM 
parameter space for n > (red crosses) and /x < (blue circles). 



gap in the plots around Aq/itiq ~ occur because it is nearly impossible to generate ruh 
as heavy as 123 — 127 GeV for such low values of trilinear couplings^. 





\i.<0 
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Figure 10: Fine-tuning measures Ahs and Aew versus Aq/too from a scan over 
mSUGRA/CMSSM parameter space for fi > (red crosses) and /i < (blue circles). 



In Fig. |11|, we plot Ahs and Aew versus tan /?. The minimal Ahs and Aew solutions 
are spread uniformly across a range of tan /? values. At very low tan /3 < 10 values, it 
is difficult to generate solutions with mh ^ 123 GeV unless mSUGRA parameters are 
extremely large, leading to high fine-tuning. 



5. Lowest fine-tuning mSUGRA benchmarks 

What is apparent from our results so far is that, after imposing LHC7 sparticle mass con- 
straints and requiring that nifi = 125 it 2 GeV on the mSUGRA/CMSSM model, the viable 
solutions are fine-tuned to at least 1% even with the less stringent EWFT measure. With 
a fine-tuning measure that knows about the high scale origin of mSUGRA parameters, the 
required fine-tuning is increased by an order of magnitude. Nonetheless, our understand- 
ing of how SUSY breaking parameters arise is extremely limited and it remains possible 
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Figure 11: Fine-tuning measures Ahs and Aew versus tan/3 from a scan over mSUGRA/CMSSM 
parameter space for > (blue circles) and < (red crosses). 



that nature may appear fine-tuned to a certain degree. With this in mind, we exhibit and 
qualitatively examine the features of the lowest Ahs and the lowest Aew solutions in the 
mSUGRA/CMSSM framework. These are listed in Table |l] as solutions HSl and EWl. 

Solution HSl has Ahs = 1100 and so requires ~ 0.1% fine-tuning. The EWFT param- 
eter Aew ~ 600, requiring ~ 0.2% fine-tuning. HSl has ttiq ~ 1500 GeV, lying at the lower 
edge of the band of solutions shown in Fig. ^. With nig ~ 1660 GeV, and rriq ~ 2000 GeV, 
this solution lies beyond the reach of LHC8 searches with up to 30 fb^^|Q, but should be 
accessible to LHC14 searches with ~ 10 — 20 fb^^ [^T|. The relatively light top squarks 
allow for g — tti decay at ~ 100%, followed by ti — tZi. Thus, gluino pair production 
will give rise to ttti+E^^^^ events at LHC and may be searchable even in the multi-jet plus 
E^^^^ channel [^]. First generation squark pair production and corresponding qg produc- 
tion will augment this rate since typically q ^ qg for first and second generation squarks. 
Production of second and third generation squarks will be suppressed by parton distribu- 
tion functions. The HSl solution has ri*~ h"^ ~ 12, so would produce too many neutralinos 
in the early universe under the standard cosmology. Late time entropy production 34 or 



neutralino decay to a lighter state, e.g. 7 + axino in extended models|3^, can bring such 
a model into accord with the measured relic abundance. The 6 — >• 57 branching frac- 
tion is somewhat below measured values, although additional flavor-violating Lagrangian 
soft terms could bring this value into accord with measurements without affecting LHC 
p henomenology. 

The solution EWl has Ahs ~ 1-5 x 10^, but Aew ~ 100 so that the latter requires 
EWFT at the 1% level. The reader may wonder whether it makes sense to talk about low 
values of Aew when Ahs is so much larger. In this connection, it may be worth allowing 
for the possibility that the mSUGRA framework may itself one day be derived from an 
underlying theory along with specific relations between seemingly unrelated mSUGRA 
parameters that lead to cancellations of the terms containing the large logarithms, as 
discussed at the end of Sec. 2^. Returning to the EWl point in the Table with nig ~ 
1600 GeV and rriq ~ 6 — 8 TeV, we see that this model is only accessible to LHC14 searches 
with ~ 50 — 100 fb^^ of integrated luminosity 1 31]. In this case, gluino pair production 
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parameter 




llvWl 


mo 


1472.0 


7926.4 


"11/2 


711.0 


594.6 




-3157.4 


5968.2 


tan p 


34.1 


29.8 




1550.8 


619.7 


m-g 


1662.5 


1589.9 


ruu^ 


2058.8 


7949.5 




2025.4 


7972.3 


mg^ 


1494.7 


7922.0 




887.8 


4547.6 




1499.8 


6197.4 




14/5.6 


61 ^^0.2 




i / oi.U 


( 4UD.D 




1 00"^ Q 


71 87 






7^fi'^ 8 
1 ouo.o 




1339.9 


7565.6 


Tfl — 


1550.1 


657.6 


777, — 


594.0 


490.4 




1547.9 


659.0 




1545.2 


638.5 


77i ~ 
^2 


591.9 


487.7 


777, ~ 
'"Zl 


308.1 


257.6 






1 9*^ 1 


5F(5 ^ S7) X 10^ 


12.3 


9.4 


2.7 


3.1 


BF{Bs fM+fi-) X 10^ 


4.4 


3.8 


a^^(Zip) (pb) 


1.4 X 10-11 


1.6 X 10-1° 


Ahs 


1105 


1.5 X 10"^ 


Aew 


582.9 


92.4 



Table 1: Input parameters and masses in GeV units for the two mSUGRA/CMSSM benchmark 
points with the lowest values of Ahs and Aew after imposing = 125 ± 2 GeV and also the 
LHC7 sparticle mass bounds. 



would be followed by gluino three-body decays to multi-jet plus multi-lepton plus E™^^^ 
final states. The final states would be rich in W and Z bosons, leading to distinctive 
signatures 1 33]. The thermally-produced neutralino abundance h? ~ 10, so again a non- 
standard cosmology as well as an extension of the spectrum is needed to bring this solution 
in accord with the measured dark matter density. 

Both HSl and EWl points will need yet other new physics to bring them in accord with 
the E821 measurement |26] of the muon magnetic moment if this discrepancy continues to 
hold up. 
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6. Concluding Remarks 



The recent discovery of a 125 GeV Higgs-like resonance at LHC has set a strong new 
constraint on super symmetric models. In addition, the lack of evidence for a SUSY 
signal at LHC now requires masses of strongly interacting sparticles in models such as 
mSUGRA/CMSSM to be above the 1 TeV scale. If LHC searches for sparticles continue 
without a new physics signal, then the little hierarchy problem - how to reconcile the Z 
and Higgs boson mass scale with the scale of SUSY breaking - will become increasingly 
acute in models such as mSUGRA.^*^ 

In this paper, we have reported on results from the calculation of two measures of fine- 
tuning in the mSUGRA/CMSSM model. The first - Ahs which includes information about 
the high scale origin of mSUGRA parameters -~ is the more stringent one. The second, Aewj 
depends only on the physical spectrum and couplings, and so is universal to all models that 
yield the same weak scale Lagrangian. Our results incorporate the latest constraints from 
LHC7 sparticle searches along with a light Higgs scalar with m/^ ~ 123 — 127 GeV. We 
find Ahs ^ 10^, or at best 0.1% fine-tuning. The more model-independent EWFT gives a 
Aew ^ 10^1 or at best 1% fine-tuning. The minimum value of ^ew tends to occur near 
the FP region which extends to large values of mo and mi/2 but which does not always 
overlap with the neutralino relic density allowed HB region. We will leave it to the reader 
to assess how much fine-tuning is too much, and also how much credence one should give 
to Ahs ill light of our ignorance of physics at or around the GUT scale^^. 

From a scan over the entire mSUGRA/CMSSM parameter space including LHC spar- 
ticle and Higgs mass constraints, we do find viable regions where EWFT is at the 1% level, 
even for gluino and squark masses well beyond LHC reach. These regions are characterized 
by mo ~ 8 TeV and ~ ib0.6mo. Since these points are spread across a wide range of 
"''1/2 values ranging up to and perhaps beyond 2 TeV, it appears that regions of parameter 
space with EWFT at the 0.5-1% levels (but with very large values of Ahs) will persist even 
after the most ambitious LHC SUSY searches are completed. 

To conclude, we remind the reader it was the realization that SUSY can solve the hig 
hierarchy problem which provided the rationale for low scale SUSY. This remains unal- 
tered by LHC and Higgs mass constraints. The underlying hope was that with sparticles 
close to the weak scale, there would be no hierarchy problem. The data seem to indicate 
that, at least in the mSUGRA framework, EWFT at the percent level is mandatory. It is 
difficult to say whether these considerations point to the failure of the mSUGRA model, 
or whether the little hierarchy is the result of an incomplete understanding of how soft 
supersymmetry breaking parameters arise. While we continue to regard models with low 
EWFT as especially interesting, it appears difficult to unilaterally discard SUSY models 
that are fine-tuned at a fraction of a percent or a part per mille, given that these provide 

^"For an illustration of how the little hierarchy problem may be solved if we go beyond mSUGRA, see 
e.g. Ref. 0. 

^^Of course, if we take the mSUGRA model to be the final high scale theory, we would no doubt take 
Ahs to be our fine-tuning measure, but the judgement to be made is whether one should treat mSUGRA 
in this manner. 



-18- 



the solution of the much more pressing big hierarchy problem. Our results provide a quan- 
titative measure for ascertaining whether or not the remaining mSUGRA/CMSSM model 
parameter space is excessively fine-tuned, and so could provide impetus for considering 
alternative SUSY models. 
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